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Several experiments have been performed on ir-TiSe 2 in order to identify whether the electronic 
structure is semimetallic or semiconducting without reaching a consensus. In this Letter, we theoret¬ 
ically study the impact of electron-hole and electron-phonon correlations on the bare semimetallic 
and semiconducting electronic structure. The resulting electron spectral functions provide a di¬ 
rect comparison of both cases and demonstrate that ir-TiSe 2 is of predominant semiconducting 
character with some spectral weight crossing the Fermi level. 


In noncorrelated materials, the semimetal or semicon¬ 
ductor character of a band-structure refers to the pres¬ 
ence or absence of bands crossing the Fermi level. In cor¬ 
related materials however, the electronic structure can 
be more subtle, showing coherent and incoherent parts 
with different spectral weights. Ab initio band-structure 
methods reducing the Hamiltonian to a single-particle 
problem result in energy-momentum relations expressed 
by ^-functions. These bare bands are not directly ac¬ 
cessible to measurements of correlated materials. Based 
on this consideration, the interpretation of the electronic 
structure should refer to the electron spectral function 
calculated via the self-energy. For lT'-TiSe 2 we have 
shown in a previous study [1] that electron-hole fluc¬ 
tuations are strong even above the transition tempera¬ 
ture, strong enough to influence the electronic structure 
around the Fermi level. 

1 T'-TiSe 2 is a layered quasi-two-dimensional (2D) com¬ 
pound with a commensurate 2x2x2 charge density 
wave (CDW) of critical temperature Tc « 190 K [2]. The 
question of the nature of the gap (positive for a semicon¬ 
ductor or negative for a semimetal) above the transition 
temperature is widely debated. In 1976, based on elec¬ 
tronic transport properties, Di Salvo et al. [5] claimed 
that lT'-TiSe 2 is a semimetal above the transition tem¬ 
perature. This assumption was confirmed in 1985 by An¬ 
derson et al. [3] performing angle-resolved photoemission 
spectroscopy (ARPES), but in 2002, Kidd et al. [1] found 
a very small indirect gap suggesting a semiconductor. In 
2007, based on optical spectroscopy, Li et al. [S] dis¬ 
agreed, claiming that their measurements clearly reveal 
that the compound is metallic in the high-temperature 
normal phase. One year later, Rasch et al. [5] chal¬ 
lenged this conclusion with ARPES, defending that the 
analysis yields undoubtedly semiconducting behavior in 
lT-TiSe 2 . 

In this Letter, we evaluate the electron-hole suscepti¬ 
bility, calculate the acoustic phonon softening driving the 
transition to a lattice deformation in the CDW phase and 
calculate the influence of Coulomb and electron-phonon 
interactions on the electron spectral function for both 
the semimetal and the semiconductor bare-band struc¬ 


ture. Starting from bare bands of either semiconducting 
or semimetallic structure we show that ARPES does not 
allow one to quantify the bare gap in the undeformed 
phase. However, there are differences between the two 
scenarios. Eirst, contrary to the semimetal, where the 
fluctuations prepare a BCS-like transition, the semicon¬ 
ductor scenario shows typical features of a Bose-Einstein 
(BE) condensation. Second, the spectral function of the 
conduction-band has a coherent peak above the Eermi 
level in the semiconductor case but not in the semimetal 
case. Third, the calculation of the phonon softening as 
a function of temperature convincingly points towards a 
bare semiconducting electronic structure. These facts al¬ 
low new insights to the semimetal versus semiconductor 
debate. The conduction-band coherent peak lays above 
the Fermi level, suggesting a semiconductor, but an inco¬ 
herent part gives a non-negligible electronic contribution 
at the Fermi level explaining the semimetallic character 
claimed in some experiments. 

Our electron-hole correlation model is based on work 
of Jerome et al. [7] discussing the ground state of an 
excitonic insulator (El) in the framework of Green func¬ 
tions. The El phase is induced by electron-hole pairs 
coupled by Coulomb interaction that condensate at a 
critical temperature, forming quasiparticles called exci- 
tons. The initial gap of the system above the transition 
temperature to the El phase determines the BCS or BE 
condensation type [HHin]. The formalism we use here is 
similar to the one developed in previous papers mm- 
Here we just mention the main steps of our calculations. 
The Hamiltonian contains a kinetic energy term com¬ 
posed by the dispersions of the topmost hole-like valence- 
band [£v(fc)] and the lowest electronlike conduction-band 
[ec(fc)] and an interaction term coupling the two bands 
via the Coulomb potential. The band dispersions are 
parabola fits to density function theory (DFT) disper¬ 
sions [12] performed with the Wien2k package [T3|. In 
the hexagonal Brillouin zone (BZ) of lT-TiSe 2 , the max¬ 
imum of the valence-band is located at the center of the 
BZ, whereas the conduction-band has its minima at its 
border. 

The extension of the El model to fluctuations above 
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the transition temperature introduces a two-particle in¬ 
teracting Green function G 2 . It is computed from the 
zeroth order in the Coulomb potential is a 

product of the Green functions for holes and electrons 
and G 2 is obtained via the Bethe-Salpeter (BS) equation 

G2 (Q,P,p','w,z) = Spp>Gi°\Q,p,w,z} 

+i^Gi°\Q,p,w,z)J2(^2(Q,P + q,p',w,z), 
q 


where w is the CDW wave vector between the centers 
of the valence and the conduction-band, and fl is the 
volume of the unit cell. In order to solve the BS equation, 
we use a local potential which is constant in the reciprocal 
space (Vb). The coordinates {Q,p,p'} are the center- 
of-mass momentum and the relative momenta. Setting 
Q = 0 is shifting the two parabola such that the centers 
match. Summed over the relative momenta p and p', 
the two-particle Green function G 2 has the properties of 
an electron-hole susceptibility [called X(Q, z), where z is 
the complex energy]. The approximation of the Coulomb 
potential by a constant Vq in the reciprocal space allows 
one to write the susceptibility X(Q,z) as a function of 
its noninteracting part X^^^Q,z) 


X(Q,z) 


X<^°\Q,z) 

l-^XW{Q^zy 


( 1 ) 


The screening length corresponding to the magnitude of 
Vq is comparable to the Thomas-Fermi screening length 
calculated for lT-TiSe 2 via the plasma frequency [S]. It 
is about one half of the nearest-neighbor distance, vali¬ 
dating the local potential approximation, as discussed in 
a previous paper [T]. The noninteracting electron-hole 
susceptibility = lim 5 _,.o w — id)] is 

numerically computed with the Cuba package |14j . The 
self-energies of electrons, and tJci and phonons Uph are 
given by 


(Tv(fc,Za) 



—X(Q + IVF(£c(fc + Q + w)) 

2tt Za + to — edk + Q + w) 


Crc(fc,Za) = 

crph{Q,Za) = lX(Q,Za) + X{Q,-Za)] ■ (2) 


Ny /B are the Fermi-Dirac and Bose-Einstein distribution 
functions, D is the local Coulomb interaction Vb/fl or 
the electron-phonon coupling g (depending on whether 
we calculate the influence of the electron-hole or electron- 
phonon correlations on the electronic structure), X is the 
spectral function of the electron-hole susceptibility X or 
of the acoustic phonon (not shown here, see Ref.[T]) and 
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FIG. 1. Electron-hole susceptibilities as a function of en¬ 
ergy for Q = 0. In blue, the susceptibilities result from a 
semimetal with a negative gap of -50 meV. In red, the calcu¬ 
lation is done for a semiconductor with a positive gap of 50 
meV. Full lines are for the exciton condensation temper¬ 
ature. Dotted lines are for > T^, the temperature of the 
phonon softening (lattice instability), a) Real and b) imagi¬ 
nary part of the noninteracting susceptibility, c) Real and d) 
imaginary part of the interacting susceptibility. 


Za is the Matsubara frequency. The renormalized phonon 
frequency is given by the condition 

- Wq - Re [crph(0,a;r. -I- i(5)] = 0, (5 0, (3) 

where tUg is the bare phonon frequency taken at the CDW 
wave vector at room temperature |15j where the self¬ 
energy correction is small. 

Here we compare these quantities starting from a semi¬ 
conducting or a semimetallic bare gap size. We set the 
chemical potential via the charge neutrality condition. 
The Coulomb potential Vg in Eq. Q is chosen such that 
X(O,0) diverges at the temperature T^. Then, we cal¬ 
culate the phonon self-energy and the phonon renormal¬ 
ization ^ in order to find a temperature for which 
tOr —>■ 0, i.e. the phonon softens, thereby inducing a lat¬ 
tice instability and the CDW phase. The parameters Vg 
and g are optimized such that the lattice instability tem¬ 
perature is equal to the experimental temperature 
of the CDW transition Tc. Finally, we calculate the elec¬ 
tron spectral functions of the valence and the conduction- 
bands. The chemical potential is set again by the charge 
neutrality condition. For bare gaps between —75 meV 
to 75 meV, the magnitude of C = Vg/^ varies between 
4.2 and 4.5 eV for semimetals and between 4.8 to 5.5 eV 
for semiconductors. The higher values for semiconduc¬ 
tors are consistent with a less-screened potential. The 
electron-phonon coupling constant g = 0.5 ± 0.1 eV/A 
for both cases and compares well with Motizuki |16j and 
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FIG. 2. (a-d) Electron spectral functions of the valence and conduction-bands for semimetal and semiconductor scenarios 
calculated via the self-energies due to electron-phonon correlations, (e-h) Same as (a-d), but with the contribution of both 
electron-hole and electron-phonon correlations to the self-energies, (i-1) Electron spectral functions corresponding to (e-h) cut 
by the Fermi-Dirac function. 


Monney et al. HI]. It turns out that as a 

consequence of Eq. ([^ being satisfied for = 0 at a 
higher temperature than the divergence temperature of 

X [Tj. 

Figure shows the electron-hole susceptibilities as a 
function of energy for the semimetal (blue curves, bare 
gap of —50 meV) and semiconductor (red curves, bare 
gap of 50 meV) cases and for two different temperatures 
at <5 = 0 . On the top, the real [Fig. [^a] and imagi¬ 
nary part [Fig. Efb] of the noninteracting susceptibility 
are presented, whereas on the bottom, the real [Fig. 
c] and the imaginary part [Fig. Sd] of the interacting 
susceptibility are shown. 

The noninteracting susceptibility X^{Q,uj) is a char¬ 
acteristic quantity calculated for CDW systems [HI US] . 
Its imaginary part [Fig. reflects the FS topology. 
The positive peak at negative energy (in blue) is char¬ 
acteristic of a semimetal band-structure and occurs at 
the energy of the negative gap. For the semiconductor 
band-structure (in red), is nonzero only for 

w larger than the gap. Its real part is the relevant quan¬ 
tity for an instability. Here the peaks are too weak to 
induce a CDW instability (about two orders of magni¬ 
tude weaker than peaks calculated for NbSe 2 , for which 
the nesting scenario is still debated [181 HO] ). 

In the interacting susceptibility X (Fig. Sc, d) a 
quasiparticle excitation occurs at very low energy. In a 
semimetal (in blue), the excitation peak grows with de¬ 
creasing temperature until diverging at the critical tem¬ 
perature (solid line) forming a long-lived quasiparti¬ 
cle, the exciton (a condensed electron-hole pair bound by 
Coulomb interaction). If the lattice was rigid (no phonon 
softening), this divergence would induce the FI phase. 


X(0, z) has approximately the form of 1 /( 2 ;^—w'^), where 
w' is the excitation energy. As in the BCS model, the 
fluctuating excitations in the electron-hole continuum are 
stable only below the condensation temperature (form¬ 
ing Cooper pairs in superconductivity and excitons here). 
For the semiconductor (in red) however, the peak already 
diverges at a non-zero energy above (e.g., at the 
dashed curves). The exciton peak lays in the gap and 
X{0,z) is well approximated by a bosonic Green func¬ 
tion of the form l/(z — oj') preparing a BE condensation 
as the temperature decreases and the peak tends towards 
a; = 0. But note that there is no exciton condensation 
here because the phonon softens at > T^. 

Figure [^presents the electron spectral functions in the 
semimetal (first column for the valence-band and sec¬ 
ond column the conduction-band) and the semiconduc¬ 
tor (third column for the valence-band and fourth column 
for the conduction-band) cases for a gap amplitude of 50 
meV. They are calculated 10 K above Tc- In the first 
row, the valence-band (VB) and conduction-band (CB) 
spectral functions A^/^{kx,ui) are shown as energy dis¬ 
persion curves (EDO), x being the direction connecting 
the center of the BZ (location of the VB) to its border 
(location of the CB). They result from the valence- and 
conduction-band self-energies due to electron-phonon in¬ 
teraction {D — g and X is the phonon spectral function). 
In the second row, the spectral functions are renormal¬ 
ized by the total self-energy i.e. the sum of the self-energy 
due to phonon interaction and due to electron-hole cor¬ 
relations where D = Vo/fl and X is the spectral function 
of the electron-hole susceptibility X. In the bottom row, 
the spectral functions displayed in false color are cut by 
the Fermi-Dirac function, showing only occupied states 
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FIG. 3. Phonon softening calculated in the semimetal (blue 
dots) and semiconductor (red dots) scenarios. These calcula¬ 
tions are compared to measurements of Holt et al. [15] (open 
black circles) and Weber et al. |21|igreen triangles). 


as accessible by ARPES. The solid lines show the bare 
bands fulfilling the charge neutrality condition and the 
dashed lines show the position of the Lorentzian fits to 
the EDC. 

Erom the comparison between the first and the sec¬ 
ond row, we conclude that the dominant effect on the 
electronic structure is due to electron-hole correlations. 
For both semimetal and semiconductor scenarios, the 
valence-band spectral function is split into two parts 
around = 0 [Figs, [^a) and[^c)]. A peak lays about 
0.15 eV above the Fermi level. As a consequence, the 
initial holelike parabola has lost some spectral weight 
around its top. This feature is qualitatively similar 
for both scenarios. For the conduction-band however, 
there is a significant difference. For the semimetal [Fig. 
[^b)], the bare conduction-band is just shifted downwards 
and crosses the Fermi level but for the semiconductor 
[Fig. id)], there is a coherent peak above the Fermi 
level (about 0.15 eV) in addition to a weaker incoherent 
peak crossing the Fermi level. If cut by the Fermi-Dirac 
function, the semiconductor coherent peak vanishes and 
the spectral functions look very similar in both scenar¬ 
ios. This similarity is a consequence of the valence-band 
spectral weight peak above the Fermi level [Figs. ga) 
and (c)]. In order to conserve the number of occupied 
states (charge neutrality), the spectral weight lost in the 
valence-band has to be compensated by electrons of the 
conduction-band below the Fermi level. Thus, electron- 
hole correlations tend to transfer some spectral weight in 
the conduct ion-band below the Fermi level independently 
of the initial bare gap (semiconductor or semimetal). 
From these considerations, it becomes clear that the com¬ 
parison with ARPES experiments does not allow one to 
decide about the nature of the initial bare gap. The cal¬ 
culated softening of the phonon as a function of temper¬ 
ature is shown in Fig. [^ together with the experimen¬ 
tal curves. For a bare semimetal (blue dots), the renor¬ 
malization is weak at high temperature and intensifies 
over a range of 10 K above the transition temperature. 
For a bare semiconductor (red dots), the renormaliza¬ 
tion is more regular. This behavior [given by Eqj^ is a 


direct consequence of the shape of the susceptibility X 
in the semimetal or semiconductor scenario illustrated in 
Fig.0 The semiconductor scenario compares much bet¬ 
ter with experiments measuring the soft phonon mode 
in ir-TiSe 2 by inelastic x-ray measurements |21j (open 
black dots) or x-ray thermal diffuse scattering [TS] (trian¬ 
gles). Note that both electron-hole and electron-phonon 
correlations cooperate in the phonon self-energy [Eq. 
to drive the transition. Though weaker, the effects of 
the electron-phonon correlations on the electron spec¬ 
tral functions tend to strengthen the effects of electron- 
hole correlations. This observation is agreement with the 
model of van Wezel et al. [22] and confirmed experimen¬ 
tally by Porer et al. |23| claiming that both the electron- 
hole and electron-phonon coupling plays a significant role 
in lT-TiSe 2 and that both effects cooperate to drive the 
CDW transition. 

Our calculations provide a first direct comparison of 
the debated semimetal versus semiconductor scenarios 
above the CDW transition temperature in lT-TiSe 2 . 
The results suggest a new way of considering the problem, 
more subtle than the conventional discussions in terms 
of bare J-functions bands. On the one hand, as sug¬ 
gested by the phonon renormalization which is directly 
proportional to the electron-hole susceptibility, this sus¬ 
ceptibility should result from a semiconductor band- 
structure. The electron spectral function is also influ¬ 
enced by the electron-hole susceptibility. The renormal¬ 
ized conduction-band shows some spectral weight cross¬ 
ing the Fermi level. These facts are not contradictory 
because this spectral weight is an incoherent part while 
the coherent peak is positioned above the Fermi level. 

This picture provides a reconciliation of the recent Let¬ 
ters 13 [6] about the nature of the gap in IT'-TiSe 2 . By 
adsorbing water on the surface, Rasch et al. |6] induced 
a band bending, the intensity of which is comparable to 
the energy position of the coherent peak above the Fermi 
level in the semiconductor case [Fig. [^d)]. Thus Rasch 
et al. deduced a semiconducting behavior from the mea¬ 
surement of the coherent peak. On the other hand, Li 
et al. 13 measured a very low carrier density using op¬ 
tical spectroscopy and concluded that the compound is 
metallic. This very low carrier density at the Fermi level 
is related to the incoherent part crossing the Fermi level 
in the semiconductor case [Fig. [^d)]. In this sense, the 
final electronic structure, taking electron-hole effects into 
account, has, at the same time, a semiconductor and a 
semimetal character. 
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